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ABSTRACT 


This  study  is  a  sensitivity  analysis  of  the  Defense  Meteorological 
Satellite  Program  (DMSP)  Block  5D  Passive  Temperature  Sounder  (SSM/T) 
to  intervening  cloud  layers.  Brightness  temperatures  for  each  of  the 
seven  SSM/T  channels  are  calculated  based  upon  two  different  climato¬ 
logical  profiles.  Intervening  precipitating  and  nonprecipitating 
cloud  layers  are  simulated  using  drop  size  distributions  based  upon 
rainfall  rate  and  Deirmendjian's  L-Model  cloud,  respectively.  Radia¬ 
tive  transfer  through  these  layers  is  approximated  using  the  discrete 
ordinate  method  for  sixteen  discrete  streams.  The  effects  of  such 
layers  on  calculated  brightness  temperatures  are  shown  to  vary  with 
cloud  position,  thickness  and  precipitation  rate.  The  importance  of 
the  variation  in  the  earth's  emissivity  between  land  and  ocean  sur¬ 
faces  is  demonstrated.  Not  onlv  do  such  variations  cause  changes  in 
brightness  temperature  values,  but  general  trends  for  individual 
channels  are  altered.  The  relative  importance  of  atmospheric,  surface 
emission  and  cloud  effects  is  broken  out  for  the  Mid-Latitude  Spring/ 
Fall  profile  over  ocean. 

The  effects  of  resulting  brightness  temperature  variations  on  the 
Air  Force  Global  Weather  Central  (AFGWC)  statistical  temperature  re¬ 
trieval  technique  and  Chahine's  relaxation  technique  are  also  presented 
The  accuracy  of  the  AFGWC  and  Chahine  temperature  retrieval  schemes  is 
demonstrated  by  comparing  the  clear  column  retrievals  to  the  actual 


profiles.  The  degradation  due  to  intervening  clouds  are  then  demon 
strated  by  presenting  the  temperature  profiles  retrieved  based  upon 
the  brightness  temperatures  calculated  for  the  cloudy  cases. 
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ABSTRACT 


This  study  is  a  sensitivity  analysis  of  the  Defense  Meteorological 
Satellite  Program  (DMSP)  Block  5D  Passive  Temperature  Sounder  (SSM/T) 
to  intervening  cloud  layers.  Brightness  temperatures  for  each  of  the 
seven  SSM/T  channels  are  calculated  based  upon  two  different  climato¬ 
logical  profiles.  Intervening  precipitating  and  nonprecipitating 
cloud  layers  are  simulated  using  drop  size  distributions  based  upon 
rainfall  rate  and  Deirmendjian's  L-Model  cloud,  respectively.  Radia¬ 
tive  transfer  through  these  layers  is  approximated  using  the  discrete 
ordinate  method  for  sixteen  discrete  streams.  The  effects  of  such 
layers  on  calculated  brightness  temperatures  are  shown  to  vary  with 
cloud  position,  thickness  and  precipitation  rate.  The  importance  of 
the  variation  in  the  earth's  emissivity  between  land  and  ocean  sur¬ 
faces  is  demonstrated.  Not  only  do  such  variations  cause  changes  in 
brightness  temperature  values,  but  general  trends  for  individual 
channels  are  altered.  The  relative  importance  of  atmospheric,  surface 
emission  and  cloud  effects  is  broken  out  for  the  Mid-Latitude  Spring/ 
Fall  profile  over  ocean. 

The  effects  of  resulting  brightness  temperature  variations  on  the 
Air  Force  Global  Weather  Central  (AFGWC)  statistical  temperature  re¬ 
trieval  technique  and  Chahine's  relaxation  technique  are  also  presented. 
The  accuracy  of  the  AFGWC  and  Chahine  temperature  retrieval  schemes  is 
demonstrated  by  comparing  the  clear  column  retrievals  to  the  actual 


profiles.  The  degradation  due  to  intervening  clouds  are  then  demon¬ 
strated  by  presenting  the  temperature  profiles  retrieved  based  upon 
the  brightness  temperatures  calculated  for  the  cloudy  cases. 
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CHAPTER  1 


INTRODUCTION 

Determination  of  atmospheric  temperature  structure  from  satellite 
sensed  atmospheric  thermal  emission  was  first  suggested  by  King  (1956) 
Kaplan  (1959)  advanced  the  concept  by  demonstrating  that  vertical  reso¬ 
lution  could  be  obtained  from  the  spectral  distribution  of  atmospheric 
emission;  observations  in  more  transparent  regions  of  the  spectrum 
sense  deeper  into  the  atmosphere.  In  order  for  atmospheric  tempera¬ 
tures  to  be  determined  by  measurements  of  thermal  emission,  the  source 
of  emission  must  be  a  relatively  abundant  gas  of  known  and  uniform 
distribution.  Carbon  dioxide,  a  minor  constituent  with  a  relative 
volume  abundance  of  .003  and  oxygen,  a  major  constituent  with  a  rela¬ 
tive  volume  abundance  of  .21,  satisfy  this  requirement.  Carbon  diox¬ 
ide  absorbs  in  an  infrared  vibrational -rotational  band  while  oxygen 
absorbs  in  a  microwave  spin-rotational  band. 

At  first  the  15  ym  COg  infrared  band  was  used  for  all  satellite 
experiments  to  determine  the  atmospheric  temperature  profiles.  The 
higher  energy  infrared  radiation  could  be  measured  by  a  much  smaller 
sensor;  an  important  consideration  for  satellite  born  systems.  From 
1963  to  1972  a  variety  of  infrared  sensors  were  flown  on  several  dif¬ 
ferent  spacecraft.  Although,  the  soundness  of  the  theory  behind  re¬ 
motely  sensing  the  temperature  profile  was  clearly  demonstrated,  one 
major  problem  continues  to  exist.  Global  atmospheric  temperature 


profile  coverage  can  not  be  obtained  due  to  the  effects  of  clouds  and 
precipitation  upon  infrared  radiation.  Several  sophisticated  techni¬ 
ques  have  been  developed  in  an  attempt  to  elliminate  these  effects 
(Smith  et.  al.,  1970;  Rodgers,  1970;  Chahine,  1970,  1974),  none  of 
which  are  valid  under  all  conditions  encountered  within  the  atmosphere. 

The  primary  advantage  of  microwave  temperature  sounders  over  in¬ 
frared  sounders  is  that  the  longer  microwaves  are  much  less  effected 
by  clouds  and  precipitation.  The  determination  of  atmospheric  temper¬ 
atures  from  satellite  measured  microwave  radiation  was  first  suggested 
by  Meeks  (1961).  Ten  years  later  the  first  Microwave  Spectrometer  ex¬ 
periment  (NEMS)  was  flown  aboard  Nimbus  5.  The  accuracy  with  which 
microwave  techniques  could  provide  vertical  temperature  profiles  was 
shown  to  be  quite  good  (Waters  et.  al.,  1975).  Furthermore,  only  0.5% 
of  all  NEMS  temperature  soundings  were  measurably  effected  by  clouds 
and  most  of  these  occured  in  the  ITCZ  (Staelin  et.  al.,  1975).  The 
major  purpose  of  the  NEMS  experiment  was  to  evaluate  passive  microwave 
techniques  for  use  on  operational  meteorological  satellites.  It  was 
not  intended  for  operational  use  itself,  and  had  only  three  channels 
which  responded  primarily  to  the  atmospheric  temperature  profile.  Of 
these  only  channel  3  (53.65GHz)  penetrated  the  atmosphere  deeply  enough 
to  be  significantly  effected  by  clouds. 

The  first  microwave  sensor  intended  for  operational  use  is  to  be 
flown  aboard  the  DMSP  Block  50  satellite  system  scheduled  to  be  launch¬ 
ed  in  1979.  This  microwave  sensor  (SSM/T)  contains  seven  channels, 
four  of  which  have  weighting  functions  peaking  near  or  below  10  km. 
These  four  channels  penetrate  deeply  enough  into  the  atmosphere  to  be 
effected  by  clouds  and  precipitation.  The  purpose  of  this  study  is  to 
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analytically  determine  what  brightness  temperature  variations  could  be 
expected  and  how  such  variations  will  impact  the  operational  tempera¬ 
ture  retrieval  scheme  to  be  used  by  the  Air  Force  Global  Weather  Cen¬ 
tral  (AFGWC) . 

Chapter  2  describes  the  characteristics  of  the  SSM/T  channels. 
Chapter  3  details  the  optical  parameters  used  and  the  cloud  and  precip¬ 
itation  models  from  which  they  were  obtained.  Chapter  4  describes  the 
solution  of  the  radiative  transfer  through  the  cloud  or  precipitation 
layer.  Chapter  5  displays  the  clear  column  calculations  for  the  two 
model  profiles  used  and  the  changes  in  brightness  temperatures  which 
result  from  clouds  and  precipitation.  Chapter  6  describes  the  AFGWC 
temperature  retrieval  scheme  and  how  the  forementioned  variations  in 
brightness  temperature  impact  this  scheme. 


CHAPTER  2 


CHARACTERISTICS  OF  SSM/T  CHANNELS 

The  SSM/T  sounder  system  operates  in  the  oxygen  absorption  band  at 
frequencies  of  50.5,  53.2,  54.35,  54.9,  58.4,  58.825  and  59.4  GHz.  The 
sensor  is  a  crosstrack  scanning  radiometer  which  acquires  data  at  32- 
second  intervals  at  7  angular  positions  separated  by  12  degrees  (see 
Figure  1.).  Channel  1  is  a  window  channel  responding  strongly  to  the 
earth's  surface  characteristics.  It  is  used  to  correct  the  other  chan¬ 
nels  for  these  background  effects.  Microwave  data  are  acquired  by 
means  of  a  mechanically  scanning,  shrouded-reflector  antenna  system. 

The  channel  parameter  design  specifications  are  presented  in  Table  1 
and  the  scan  parameters  are  illustrated  in  Table  2.  The  seven  SSM/T 
channels  have  been  selected  such  that  the  weighting  functions  peak  at 
different  altitudes  within  the  atmosphere  providing  the  necessary  ver¬ 
tical  resolution  for  temperature  profile  retrieval.  The  weighting 
functions  for  each  channel  are  displayed  in  Figure  2.  These  weighting 
functions  result  from  consideration  of  atmospheric  absorption  due  to 
water  vapor  and  molecular  oxygen  with  transmissivity  corrections  for 
antenna  gain  characteristics  included.  The  antenna  gains  employed  for 
this  purpose  are  presented  in  Table  3  and  resulted  from  an  interpola¬ 
tion  of  percent  gain  vs  deviation  from  boresight  to  the  desired  bore 
sight  deviation  angles.  This  was  accomplished  for  each  of  the  seven 
bore  sight  directions.  The  precent  gain  for  each  deviation  angle  was 
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then  averaged  over  the  seven  bore  sight  directions  to  obtain  mean  val 
ues  for  each  channel.. 


Cold  Sky 

Calibration 

Position 


Hot  Load 

Calibration 

Position 


Figure  1.  SSM/T  Scan  Geometry. 

Source:  Riogone  and  Stogryn  (1977) 

Table  1.  Channel  Parameter  Design  Specifications 


Channel  Polarization 

Frequency 

(GHz) 

Bandwidth 

(MHz) 

NETD 

(°K) 

1  -i 

50.5 

400 

0.6 

2 

Principally 

53.2 

400 

0.4 

3 

Horizontal 

54.35 

400 

0.4 

4  J 

54.9 

400 

0.4 

51 

i 

Orthogonal 

58.4 

115 

0.5 

6 

to  Channels 

1-4 

58.825 

400 

0.4 

7  J 

59.4 

250 

0.4 

Source:  Riogone  and  Stogryn  (1977) 
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Table  2.  Scan  Parameters 


Scan  Type  . 

Cross-Track  Positions  . 

Calibration  Positions  . 

Instantaneous  Field  of  View  .  .  . 

Total  Cross-Track  Scan  . 

Total  Scan  Period  . 

Integration  Time  (Crosw*  £k 
and  Calibration  Positions)  .  .  . 

Sync  Mode .  . 


Cross-Track  Nadir 
7 

2-Cold  Space  and  300°K 
14° 

36° 

32  Seconds 

2.7  Seconds 

Auto  Mode  of  On-Sync 


Source:  Riogone  and  Stogryn  (1977) 
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Table  3.  Percent  Gains  Averaged  Over  Bore  Sight  Angles 


Channel  Frequencies  (GHz)  j 


0° 

50.50 

53.20 

54.35 

54.90 

58.40 

58.82 

59.40 

0. 

99.357 

99.675 

99.675 

99.675 

100.000 

99.727 

99.350 

1. 

97.583 

97.901 

97.901 

97.901 

97.360 

96.987 

96.473 

2. 

93.819 

93.877 

94.217 

94.379 

91.545 

91.577 

91.622 

3. 

88.301 

88.095 

88.772 

89.096 

83.173 

83.012 

82.790 

4. 

80.706 

79.643 

80.360 

80.702 

72.161 

71.892 

71.519 

5. 

70.741 

68.905 

69.138 

69.249 

59.177 

58.694 

58.028 

6. 

59.742 

56.612 

56.192 

55.991 

46.061 

45.185 

43.976 

7. 

48.659 

43.934 

43.228 

42.890 

34.283 

33.306 

31.956 

8. 

37.555 

32.230 

31.664 

31.393 

24.283 

23.550 

22.536 

9. 

27.325 

22.201 

21.861 

21.698 

15.793 

15.246 

14.490 

10'. 

18.778 

13.867 

13.930 

13.961 

8.990 

9.082 

9.209 

11. 

12.523 

6.893 

6.926 

6.942 

5.466 

5.266 

•  4.990 

12. 

7.698 

4.729 

4.694 

4.678 

2.764 

2.936 

3.173 

13. 

5.189 

3.129 

2.899 

2.789 

1.691 

1.613 

1.505 

14. 

3.539 

1.371 

1.495 

1.316 

.899 

.846 

.773 

15. 

2.212 

1.106 

.725 

.543 

.441 

.421 

.393 

16. 

1.353 

.601 

.335 

.208 

.311 

.288 

.255 

17. 

.834 

.282 

.220 

.190 

.229 

.216 

.199 

18. 

.464 

.247 

.230 

.223 

.166 

.163 

.158 

19. 

.267 

.249 

.225 

.213 

.130 

.146 

.167 

20. 

.194 

.230 

.197 

.182 

.096 

.125 

.166 

21. 

.157 

.192 

.164 

.150 

.064 

.096 

.141 

22. 

.123 

.163 

.140 

.128 

.048 

.077 

.118 

23. 

.093 

.137 

.118 

.108 

.051 

.070 

.096 

24. 

.070 

.113 

.098 

.091 

.049 

.062 

.080 

25. 

.048 

.093 

.085 

.081 

.043 

.055 

.072 

26. 

.034 

.073 

.072 

.071 

.037 

.049 

.065 

27. 

.038 

.058 

.061 

.063 

.033 

.043 

.058 

28. 

.042 

.048 

.052 

.054 

.030 

.039 

.051 

29. 

.039 

.039 

.044 

.046 

.028 

.034 

.042 

Table  3.  (Continued) 


Channel  Frequencies  (GHz) 


0° 

50.50 

53.20 

54.35 

54.90 

58.40 

58.82 

59.40 

30. 

.032 

.029 

.035 

.038 

.025 

.028 

.032 

31. 

.025 

.025 

.029 

.030 

.024 

.025 

.025 

32. 

.019 

.022 

.025 

.026 

.023 

.021 

.017 

33. 

.013 

.024 

.027 

.029 

.022 

.018 

.013 

34. 

.011 

.026 

.031 

.033 

.022 

.018 

.012 

35. 

.010 

.027 

.033 

.035 

.021 

.017 

.012 

36. 

.011 

.027 

.030 

.032 

.020 

.017 

.013 

37. 

.012 

.025 

.027 

.028 

.020 

.017 

.013 

38. 

.014 

.024 

.024 

.024 

.019 

.017 

.013 

39. 

.016 

.021 

.020 

.019 

.019 

.016 

.013 

40. 

.019 

.019 

.016 

.015 

.018 

.016 

.015 

41. 

.021 

.018 

.014 

.012 

.017 

.016 

.016 

42. 

.021 

.017 

.012 

.010 

.016 

.016 

.017 

43. 

.021 

.018 

.011 

.008 

.015 

.016 

.018 

44. 

.021 

.018 

.010 

.006 

.014 

.016 

.019 

45. 

.020 

.017 

.009 

.005 

.013 

.016 

.020 

46. 

.020 

.015 

.008 

.005 

.013 

.016 

.021 

47. 

.019 

.014 

.008 

.006 

.012 

.016 

.021 

48. 

.018 

.013 

.009 

.007 

.012 

.016 

.022 

49. 

.017 

.012 

.009 

.008 

.012 

.016 

.022 

50. 

.016 

.012 

.010 

.009 

.011 

.016 

.022 

51. 

.016 

.011 

.011 

.011 

.011 

.015 

.021 

52. 

.015 

.010 

.011 

.011 

.010 

.014 

.021 

53. 

.015 

.009 

.011 

.012 

.009 

.013 

.019 

54. 

.014 

.009 

.011 

.012 

.009 

.012 

.018 

55. 

.014 

.008 

.011 

.012 

.008 

.011 

.016 

56. 

.013 

.008 

.011 

.012 

.007 

.010 

.014 

57. 

.013 

.007 

.010 

.012 

.007 

.009 

.012 

58. 

.012 

.007 

.010 

.011 

.006 

.008 

.011 

59. 

.011 

.007 

.009 

.010 

.006 

.007 

.009 

CHAPTER  3 


OPTICAL  PARAMETERS  FOR  SSM/T  CHANNELS 


3.1  Mie  Theory 

The  application  of  Mie  scattering  theory  enables  calculation  of 
optical  parameters  for  assumed  spherical  particles;  given  the  complex 
index  of  refraction,  wave  length  of  radiation  to  be  scattered  and  par¬ 
ticle  size  density  distribution.  The  Mie  program  used  for  these  calcu¬ 
lations  is  a  modification  for  Univac  computers  of  the  program  by  Liou 
and  Hansen  (1971).  The  output  of  the  Mie  computations  include  the  co¬ 
efficients  of  the  Legendre  polynomial  approximation  to  the  phase  func¬ 
tion,  the  extinction  coefficient  (Km“^)  and  the  single  scattering  albedo. 
Note  that  the  extinction  coefficient  and  single  scattering  albedo  deter¬ 
mined  here  include  scattering  and  absorption  by  the  particles  as  the 
only  mechanisms  for  extinction.  These  values  are  modified  to  include 
absorption  by  the  gaseous  atmosphere  within  the  cloud  prior  to  their 
use  in  the  solution  of  the  radiative  transfer  equation. 

3.2  Indices  of  Refraction 

The  complex  indices  of  refraction  for  pure  liquid  water  at  temper¬ 
atures  of  +10°  and  -10°  C  have  been  determined  by  Savage  (1976)  for 
several  different  frequencies.  These  values  (Table  4)  are  based  on  a 
least  squares  fit  to  the  data  of  Saxton  and  Lane  (1962)  given  by 
Hollinger  (1973).  The  complex  indices  of  refraction  for  the  SSM/T 
channels  (Table  5)  were  obtained  from  Table  4  by  linear  interpolation 
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Table  4.  Refractive  Index  m=m  -im.  of  Pure  Liquid  Water 

r  i 


Frequency 

T=10°C 

T=0°C 

T=-10°C 

(Ghz) 

mr  m. 

mr  mi 

mr 

19.35 

5.87 

2.96 

5.20 

2.94 

4.64 

2.81 

37.00 

4.32 

2.60 

3.83 

2.37 

3.46 

2.12 

50.30 

3.76 

2.29 

3.37 

2.03 

3.08 

1.78 

89.50 

3.00 

1.67 

2.76 

1.42 

2.60 

1.20 

100.00 

2.89 

1.56 

2.68 

1.31 

2.54 

1.10 

118.00 

2.75 

1.39 

2.58 

1.15 

2.46 

0.96 

130.00 

2.68 

1.30 

2.53 

1.07 

2.43 

0.89 

183.00 

2.49 

1.00 

2.39 

0.81 

2.33 

0.66 

231 .00 

2.40 

0.82 

2.33 

0.66 

2.29 

0.53 

After  Savage  (1976) 


Table  5.  Refractive  Index  of  Pure  Liquid  Water  (SSM/T) 


Channel  Frequencies  (GHz) 


50.50 

53.20 

54.35 

54.90 

58.40 

58.82 

59.40 

3.37 

3.32 

3.31 

3.30 

3.24 

3.24 

3.23 

2.03 

1.98 

1.97 

1.96 

1.90 

1.90 

1.89 

in  frequency.  For  the  purpose  of  this  study  the  values  for  0°C  are 
used  throughout.  It  is  important  to  note,  however,  that  after  the 
Block  50  satellite  flys,  greater  care  may  be  necessary  to  insure  that 
the  indices  of  refraction  are  compatible  with  atmospheric  temperatures 
at  cloud  levels  if  observed  brightness  temperatures  are  to  be  matched 
by  theoretical  calculations. 


3.3  Drop  Size  Distributions 

Once  the  absorption  and  scattering  efficiences  for  single  parti¬ 
cles  have  been  determined  by  Mie  theory  in  terms  of  effective  cross 
sectional  area,  the  efficiencies  for  a  distribution  of  particles  can  be 
determined  by  integrating  over  the  distribution.  This  area  per  volume 
relationship  then  has  dimensions  of  inverse  length  and  is  described  as 
the  absorption  (or  scattering,  or  extinction)  coefficient.  The  drop 
size  distributions  used  for  this  study  are  necessarily  numerical  models. 
Two  models  have  been  chosen  to  represent  nonprecipitating  and  precipi¬ 
tating  clouds,  respectively.  Deirmendjian's  L-Model  cloud  (Deirmendjian, 
1969)  drop  size  distribution  is  used  for  the  nonprecipitating  case,  that  is: 

N(r)  =  araexp(-br5),  (3-1) 

where  r  is  the  drop  radius  in  mm.  By  varying  the  constants  this  model 
can  be  made  to  represent  a  variety  of  loud  types.  For  this  study, 
a=4.9757X107,  a=2,  b=15.1186,  and  6=.5.  For  the  precipitating  case,  a 
theoretical  drop  size  distribution  based  on  rainfall  rate  is  used 
(Marshall  and  Palmer,  1948).  The  distribution  is  expressed  by: 

N(r)  a  0.08  exp(-2r6)  (3-2) 

where  <S=41R“‘21  and  R  denotes  rainfall  rate  (mm/hr.).  This  exponential 
behavior  of  drop  size  distribution  has  been  experimentally  verified  by 
Gunn  and  Marshall  (1958)  and  Sekhon  and  Srivastaval  (1970). 

3.4  Parameter  Values 

The  Legendre  polynomial  phase  function  coefficients,  extinction 
coefficients  and  single  scattering  albedos  calculated  for  each  SSM/T 
channel  from  Mie  theory  are  tabulated  in  Tables  6-16.  Since  the 


extinction  coefficient  input  to  the  radiative  transfer  program  must 
be  in  units  of  (KnT^)  the  final  values  ouput  by  the  Mie  program  have 
been  multiplied  by  the  proper  constant  to  obtain  these  units. 
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CHAPTER  4 


MICROWAVE  RADIATIVE  TRANSFER 


4.1  Clear  Column  Brightness  Calculations 

For  the  purpose  of  calculating  the  brightness  temperature  of  a 
given  microwave  channel  at  the  top  of  the  earth's  atmosphere,  it  is 
assumed  that: 

a)  the  atmosphere  is  in  local  thermodynamic  equilibrium, 

b)  the  atmosphere  is  plane  parallel  and  horizontally  homogenious, 

c)  the  atmosphere  is  axially  symmetric  about  nadir-zenith,  and 

d)  there  is  no  source  of  illumination  exterior  to  the  earth-atmos¬ 
phere  system. 

The  basic  equation  of  transfer  for  an  absorbing  medium  can  be 
written : 


dl 

V  —  T  1 

k  pds  v  v 
v 


(4-1) 


Assuming  transfer  along  a  path  perpendicular  to  the  surface  of  the 

earth,  dz  can  be  substituted  for  ds.  The  formal  solution  of  Eq.  (4-1) 

is :  2  z 

-Jk^pdz'  2  /kvpdz' 

I  (z)  =  I  (0)e  0  V  +  /  J  (z')e  0  k  pdz*  (4-2) 

v  v  0  v  v 

where  ky  is  the  coefficient  of  gaseous  absorption  and p  is  the  density  of 

absorbing  gases.  Defining  the  transmissivity  as: 

z 

T  (z‘  ,z)  =  exp(-J  k  pdz), 

V  £  i  V 


(4-3) 
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enables  equation  (4-2)  to  be  written: 

I  (z)  =  I  ( 0 ) t  (0,z)  +  /  J  (z')dT  (z',z). 

v  v  v  z‘=0  v  v 

Therefore,  the  radiant  energy  reaching  the  top  of  the  atmosphere  (z=®) 
is  given  by: 

Iv  =  Iv(0)tv(0,®)  +  J  Ov(z)dxv(z' ,®)  (4-4) 

Under  assumption  (a)  the  radiant  energy  emitted  at  any  point  in 
the  atmosphere  is  a  function  of  the  temperature  and  frequency  and  is 
given  by  Planck's  Function.  In  the  frequency  domain  the  Planck  Func¬ 
tion  is: 

Bv(T)  =  2hv3/[c2(ehv/KT-l)]  (4-5) 

where  v  is  frequency,  h  is  Planck's  constant,  c  is  the  speed  of  light 
and  K  is  Boltzman's  constant.  In  the  microwave  region  hv/KT«l, 
therefore. 


exp(hv/KT)  =  1+hv/KT. 


Then  equation  (4-5)  reduces  to 


B  (T)  =  (^A  T  =  CT. 


(4-6) 


(4-7) 


In  other  words  for  a  given  frequency  the  Planck  Function  is  simply  a 
constant  times  the  absolute  temperature.  Moreover,  the  equivalent 
brightness  temperature  for  upwelling  intensity  I  may  then  be  defined 
as 


I 

v 


CTb(v). 


(4-8) 
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The  radiant  energy  emitted  at  a  given  level  is  isotropic  and  makes 
two  contributions  to  the  radiant  energy  reaching  the  top  of  the  atmos¬ 
phere  : 

a)  energy  emitted  in  the  upward  direction  travels  directly  to  the 
top  of  the  atmosphere  while  suffering  a  Beer's  Law  attenuation  along 
its  path  length  due  to  absorption  by  atmospheric  gases, 

b)  energy  emitted  in  the  downward  direction  is  attenuated  by  ab¬ 
sorption  along  its  path  length  to  the  surface,  where  a  portion  of  it  is 
reflected  back  upward  and  then  attenuated  along  its  path  length  to  the 
top  of  the  atmosphere.  Mathematically,  the  contribution  (a)  reaching 
the  top  of  the  atmosphere  from  height,  z,  can  be  expressed  as: 

T |J(v,z)  *  T(z)tv(z,»).  (4-9) 

The  contribution  (b)  reaching  the  top  of  the  atmosphere  is  given  by: 

TjJ(v.z)  =  T(z)(l-evhv(0,-)Tv(0,z).  (4-10) 

Normally,  transmissivities  are  calculated  as: 

t  (z ,<*>)  »  exp (-/  k  pdz' ) ,  (4-11) 

z 

representing  the  portion  of  the  upwelling  radiance  at  level  z  which 
reaches  the  top  of  the  atmosphere.  It  can  be  easily  shown  from  Eq. 
(4-3)  that: 

tv(0,z)  *  tv(0,»)/tv(z,®)  (4-12) 

Similarly, 

Tv(z\z)  =  tv(z' ,«)/tv(z,») 


4 


(4-13) 
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for  any  layer  defined  by  z'  and  z(z'<z). 

To  obtain  the  total  atmospheric  contribution  to  radiance  reaching 
the  top  of  the  atmosphere  in  any  channel,  simply  integrate  equations 
(4-9)  and  (4-10)  and  add  them  together: 

Tj(v)  =  T°(v)  +  Tj(v)  =  /  T(z)dx  (z,») 

°  D  0  z=0 


+  (1-e  )t  (0,»)  /  T(z)dx  (0,z) 


This  can  be  approximated  numerically  as: 


Tb(v)  =  ETiAi.Tv(z,®)  +  (1-ev)tv(0,»)  lTiA1-Tv(0,z) 


(4-14) 


(4-15) 


The  surface  of  the  earth  emits  in  the  microwave  range  as  a  grey 
body.  Therefore,  the  background  emission  of  the  earth's  surface  atten¬ 
uated  to  the  top  of  the  atmosphere  by  Beer's  Law  is  given  by: 


TB(v)  =  evVv<0*-> 


Combining  equations  (4-15)  and  (4-16), 


(4-16) 


Tg(v)  =  [evT$  +  zTi(l-ev)AiTv(0,z)]Tv(0,®)  +  TT^t^Z,-*)  .  (4-17) 


Now  note  that: 


^iTv(0,z)  =  Tv(0»z.j  )-xv(0,z1.+1 ) 


Tv(0,«)  Tv(0,®) 


•m>  VTi+1*" 


T  (Z, 

V  1 


Tv(0,«)[xv(zi+1 ,«)-Tv(Zi ,«)] 


TvUi’“,Tv(2i+l*“) 
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’oo)Tv(zi+r°,)J 


A1Tv^z**) 


Therefore  equation  (4-17)  can  be  written: 


¥*> *  *vV,<0-> 


n-1  t  (0,«) 

+  i-i  C(,*I'>)  *  ,]Vi\,<2-"> 


(4-18) 


where  n  is  the  number  of  levels  used  to  represent  the  atmosphere.  This 
is  simply  the  approximation  to: 


tB(v)  =  evT$Tv(0,«)  +  /  [(l-ev)tv(0,z)Z+l]T(z)dtv(z.-)  (4-19) 


Equation  (4-19)  is  the  absorption/emission  solution,  equation  (4-4),  to 
the  radiative  transfer  equation  where. 


I  (0)  -  C  e  T5  (4-20) 

V  V  S  '  9 

and 

Jv(z)  =  C  T(z)  [(l-ev)xv(0,z)2  +  1]  (4-21) 

This  is  a  good  approximation  for  a  clear  atmosphere  where  scattering 
in  the  microwave  region  is  very  small.  Note  that  equation  (4-19) 
yields  what  is  commonly  referred  to  as  antenna  temperature  if  the 
transmissivities  used  have  been  corrected  to  include  antenna  gain 
characteristics.  With  unaltered  transmissivities  equation  (4-19) 
yields  the  clear  column  brightness  temperature. 


i 


4.2  Brightness  Calculations  for  Cloudy  Atmospheres 

The  basic  equation  of  transfer  for  a  plane  parallel  cloud  layer 
in  local  thermodynamic  equilibrium  can  be  written  in  the  form: 


udIv(t,y)/dT  =  I v( t ,p )-Cov/2  j  Pv(u,u‘ )Iv(T,ii' )du' 


-  (1-5>v)Bv[T(t)]  (4-22) 

where  1^  represents  the  monochromatic  radiance,  u  the  cosine  of  the 
emergent  angle  with  respect  to  the  zenith,  t  the  monochromatic  opti¬ 
cal  depth,  Py  the  normalized  axially  symmetrical  scattering  phase 
function  and  the  single  scattering  albedo.  Bv[T(x)]  is  the  Planck 
Function  as  defined  by  equation  (4-7).  The  normalized  phase  function 
can  be  expressed  as  a  Legendre  polynomial  of  finite  terms.  By  approx¬ 
imating  the  integration  in  equation  (4-22)  utilizing  Gauss'  quadrature 
formula,  a  set  of  first  order  inhomogeneous  differential  equations 
can  be  derived.  The  homogeneous  and  particular  solutions  of  the  dif¬ 
ferential  equations  as  outlined  by  Chandrasekhar  (1950),  enables  the 
complete  solution  of  the  scattered  intensity  from  an  isothermal  cloud 
(temperature,  T  )  for  a  given  discrete  stream  to  be  expressed  as: 

iv<T’“i> =  W“i>exP<-V>t|VTc>  l4-23> 

m 

where  m  is  the  number  of  discrete  streams  employed,  0m  and  km  represent 
eigen  functions  and  eigen  values  for  the  differential  equations  and 
are  a  set  of  constants  of  proportionality  which  can  be  determined  from 
the  radiation  boundary  conditions.  Noting  that  the  brightness  tempera¬ 
ture  in  the  microwave  range  is  simply  the  radiant  intensity  derived  by 
a  constant,  equation  (4-22)  can  also  be  written  in  terms  of  brightness 
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temperature.  The  associated  solution  is  then  given  by: 

Tb(v,t;u.)  =  ILm0m(^)exP(-kmT)+Tc  (4"24) 

m 

where  the  constants  of  proportionality,  Lm,  are  different  from  those  in 
equation  (4-23),  having  been  obtained  using  brightness  boundary  condi¬ 
tions.  T  is  the  cloud  layer  temperature  and  is  independent  of  fre- 
quency,  v. 

At  the  cloud  top,  the  downward  brightness  is  equal  to  the  bright¬ 
ness  contributions  from  every  point  in  the  atmosphere  above  the  cloud 
top  attenuated  by  Beer's  Law  to  the  cloud  top.  This  can  be  expressed 
as : 

Vv»2t;'wi>  =  7  T(z)  dTv(z,zt;"ui)  (4~25) 

z=zt 

where  z^  is  the  height  at  the  cloud  top  and  the  negative  sign  on 
simply  indicates  downward  transfer.  Note  that: 

xv(zt,z;y.)  «  Tv(z,zt;-yi)  =  expH^-t)/^]  (4-26) 

and 

t v(z1,z2;y.)  *  Tv(z1,»;y.)/Tv(z2,«;y.)  (4-27) 

where  z2>z-| 

Within  the  cloud  layer,  where  scattering  occurs,  continuity  of 
brightness  values  from  all  directions  is  required.  Thus: 

Tfi ( v»Z£ »y.j )  =  ^B^V*ZA+1  *wi ^  i_l  ,2 , . . .  ,N-1  (4-28) 

where  N  is  the  total  number  of  sublayers  within  the  cloud. 
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At  the  lower  boundary  of  the  cloud,  three  brightness  contribu¬ 
tions  are  immediately  apparent: 

a)  the  surface  contribution, 

b)  the  direct  atmospheric  contribution  from  below  the  cloud  and 

c)  the  reflected  atmospheric  contribution  from  below  the  cloud. 

A  fourth,  and  perhaps  less  obvious,  brightness  contribution  at  the 
lower  boundary  must  be  considered.  Since  the  reflection  of  downward 
brightness  by  the  earth's  surface  is  significant  for  microwave  radia¬ 
tion,  the  emergent  brightness  at  the  cloud  bottom  will  contribute  to 
the  lower  boundary  condition.  That  is  to  say,  the  solution  to  the 
radiative  transfer  through  the  cloud  effects  the  boundary  conditions 
used  to  obtain  the  solution.  This  suggests  that  an  iterative  approach 
to  the  correct  solution  is  required.  Practically,  this  is  accomplished 
by  assuming  initially  that  the  top/down  throughput  of  the  cloud  equals 
one.  The  emergent  intensity  at  the  cloud  bottom  obtained  by  this 
assumption  can  then  be  attenuated  by  two  trips  through  the  atmosphere 
between,  the  earth's  surface  and  the  cloud  base  and  by  the  reflection 

at  the  earth's  surface.  This  value  is  added  to  the  lower  boundary 
condition  and  the  solution  to  the  radiative  transfer  equation  then 
produces  a  new  value  for  the  emergent  brightness  at  the  cloud  bottom. 
From  this  solution  a  new  top/down  throughput  can  be  calculated  and  the 
process  is  repeated  until  the  new  throughput  varies  from  the  old  by 
less  than  one  tenth.  In  all  cases  examined  by  this  study  the  itera¬ 
tion  halted  after  only  two  steps.  The  lower  boundary  condition  can 
then  be  expressed  as: 
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Tb( v,zb ;vi )  =  evTsTv(0,zbiyi^  +  zlQ  T(z)dTv(z,zb;Ui) 


+  (l_ev)Tv(0,zfc);u1.)[  /  T(z)c1tv(0,z;-p1.) 


+  TB^v,zt;“ui^Tp^yi^Tv^0*zb’vi^ 


(4-29) 


where  Tp( u-f )  is  defined  to  be  the  top/down  throughput  of  the  cloud  for 
the  stream  defined  by  and  zb  is  the  height  of  the  cloud  base. 


Tp(yi)  =  TB^v’zb’"ui^TB^v’zt’"ui^ 


(4-30) 


The  solution  of  the  microwave  transfer  equation  given  by  equation 
(4-24)  is  applicable  only  to  isothermal  and  homogeneous  cloud  layers. 

To  apply  the  transfer  solution  to  thick  clouds  it  is  necessary  to  di¬ 
vide  the  cloud  into  sublayers  for  which  the  isothermal  and  homogeneous 
assumptions  are  reasonable.  By  matching  the  brightness  continuity 
equation  for  brightness  in  each  sublayer,  a  set  of  linear  equations 
with  unknown  coefficients  can  be  obtained  which  can  be  solved  by  stand¬ 
ard  matrix  inversion  techniques.  Similar  procedures  have  been  used  by 
Liou  (1975)  to  evaluate  the  transfer  of  solar  radiation  in  inhomogen¬ 
eous  atmospheres  and  by  Feddes  and  Liou  (1977)  to  investigate  the 
transfer  of  spectral  infrared  radiation  in  cloudy  atmospheres. 


CHAPTER  5 


SENSITIVITY  TO  CLOUD  LAYERS 

Sections  5.1  through  5.4  of  this  chapter  are  descriptive  in  nature. 
The  profiles  and  cloud  models  used  are  discussed  in  section  5.1  Sec¬ 
tions  5.2  through  5.5  simply  report  the  results  obtained  with  no 
attempt  at  explanation.  Section  5.5  deals  with  the  interpretation  of 
the  results  and  the  physical  mechanisms  which  produced  them. 

5.1  Profiles  and  Cloud  Models  Used 

Two  climatological  profiles  were  used  for  this  sensitivity  study. 
The  first  profile  was  derived  from  the  Northern  Hemispheric  Mid-Lati¬ 
tude  Spring/Fall  climatology  contained  in  the  U.S.  Standard  Atmospheric 
Supplements,  1966.  The  second  profile  was  derived  from  the  climatology 
for  30°  North  Latitude,  July,  contained  in  the  same  source.  For  both 
profiles  temperature  and  height  values  were  interpolated  to  40  standard 
pressure  levels  used  for  this  study.  Mixing  ratio  values  for  levels 
below  250  mb  were  obtained  from  a  Skew-T  Diagram  and  hence  are  satura¬ 
tion  values.  Above  250  mb  a  constant  mixing  ratio  of  3  parts  per 
million  was  used.  The  resulting  Mid-Latitude  Spring/Fall  profile  is 
given  in  Table  17,  while  the  30°  N.  Latitude,  July,  profile  is  contain¬ 
ed  in  Table  18. 

Cloud  models  of  1,  2,  3,  4  and  5  Km  thicknesses  were  investigated, 
employing  the  Marshall -Palmer  drop  size  distribution  to  simulate  rain¬ 
fall  rates  of  1,  2,  3,  4,  5,  10,  15,  20,  25  and  30  mm/hr.  In  addition. 


Table  17*  N.  Hemisphere  Mid-Latitude  Sprinq/Fall  Profile 


Pres(mb) 

Height(Km) 

Temp(°K) 

Vapor (gm/m^) 

1000.00 

.no 

287.42 

.1240E+2 

950.00 

.540 

284 . 64 

.  1049E+2 

920.00 

.810 

282.90 

. 9248E+1 

899.54 

1.000 

281.70 

.8625E+1 

850.00 

1.460 

278.67 

. 71 38E+1 

794.81 

2.000 

275.12 

. 5682E+1 

780.00 

2.150 

274.13 

•5298E+1 

700.89 

3.000 

268.64 

- 3530E+1 

700.00 

3.010 

268.57 

•3511E+1 

670.00 

3.360 

266.32 

.3094E+1 

620.00 

3.960 

262.43 

. 2254E+1 

616.78 

4.000 

262.17 

.2213E+1 

570.00 

4.600 

258.24 

-1625E+1 

540.60 

5.000 

255.66 

-1349E+1 

500.00 

5.580 

251 .85 

. 9300E+0 

475.00 

5.960 

249.47 

-8078E+0 

472.41 

6.000 

249.21 

•7920E+0 

430.00 

6.680 

244.77 

•5391E+0 

400.00 

7.190 

241.44 

- 3875E+0 

356.57 

8.000 

236.19 

-2374E+0 

350.00 

8.130 

235.34 

-2154E+0 

300.00 

9.180 

228.54 

.1061E+0 

265.03 

10.000 

223.23 

. 6446E-1 

250.00 

10.380 

220.73 

•4735E-1 

200.00 

11.810 

216.65 

.9650E-3 

194.29 

12.000 

216.65 

.9374E-3 

150.00 

13.650 

216.65 

.7237E-3 

135.00 

14.310 

216.65 

•6513E-3 

115.00 

15.330 

216.65 

•5549E-3 

100.00 

16.220 

216.65 

•4825E-3 

85.00 

17.260 

2)6.65 

.4101E-3 

70.00 

18.500 

216.65 

•3377E-3 

60.00 

19.480 

216.65 

. 2895E-3 

50.00 

20.640 

217.26 

•2406E-3 

30.00 

23.940 

220.51 

•1422E-3 

25.00 

25.130 

221.69 

.1179E-3 

20.00 

26.590 

223.14 

•9369E-4 

15.00 

28.500 

225.02 

-6968E-4 

10.00 

31.210 

227.18 

•4601E-4 

7.00 

33.630 

224.11 

•3265E-4 

5.00 

35.980 

239.14 

•2186E-4 

4.00 

37.580 

243.94 

•1714E-4 

3.00 

39.680 

249.62 

•1256E-4 

2.00 

42.730 

258.02 

.8102E-5 

1.50 

44.950 

264.05 

-5938E-5 

1.00 

48.180 

270.52 

-3864E-5 

.50 

53.760 

267.65 

•1953E-5 

.20 

60.880 

253.96 

.8232E-6 

.10 

65.960 

235.47 

•4439E-6 

I 

I 


Table  18.  30°  N.  Latitude  July  Profile 


37 


1 

1 

1 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 


Pres(mb) 

Height(Km) 

Temp(°K) 

H^O  Vapor(gm/m^) 

1000.00 

.120 

303.51 

. 3107E+2 

950.00 

.570 

299.44 

.2411E+2 

920.00 

.850 

297.12 

. 2182E+2 

904.78 

1.000 

296.12 

.2071E+2 

850.00 

1.540 

292.36 

. 1 676E+2 

805.30 

2.00 

289.61 

.  1424E+2 

780.00 

2.270 

287.99 

. 1 283E+2 

714.56 

3.000 

283.73 

.9815E+1 

700.00 

3.170 

282.73 

. 91 58E+1 

670.00 

3.540 

280.60 

.8115E+1 

633.18 

4.000 

277.92 

. 6794E+1 

620.00 

4.170 

276.92 

•6329E+1 

570.00 

4.850 

273.05 

•4914E+1 

559.45 

5.000 

272.21 

-4655E+1 

500.00 

5.890 

267.16 

.3235E+1 

493.03 

6.000 

266.43 

. 3049E+1 

475.00 

6.290 

264.51 

. 2573E+1 

430.00 

7.050 

259.10 

. 1 738E+1 

400.00 

7.600 

255.25 

. 1 257E+1 

379.05 

8.000 

252.45 

. 1024E+1 

350.00 

8.590 

248.31 

.7105E+0 

300.00 

9.700 

240.53 

. 3607E+0 

287.37 

10.000 

238.44 

•3062E+0 

250.00 

10.960 

231.66 

.1514E+0 

214.13 

12.000 

224.48 

.7763E-1 

200.00 

12.450 

221.31 

. 5005E-1 

150.00 

14.260 

208.69 

. 51 59E-2 

135.00 

14.910 

204.20 

.6911E-3 

115.00 

15.870 

203.20 

. 5916E-3 

100.00 

16.710 

204.56 

•5110E-3 

85.00 

17.690 

206.71 

.4298E-3 

70.00 

18.889 

209.31 

.3496E-3 

60.00 

19.840 

211.40 

.2967E-3 

50.00 

20.980 

213.89 

.2444E-3 

30.00 

24.260 

220.41 

.1423E-3 

25.00 

25.450 

222.78 

•1173E-3 

20.00 

26.930 

225.71 

.9262E-4 

15.00 

28.860 

229.54 

.6831E-4 

10.00 

31.650 

235.10 

•4446E-4 

7.00 

34.160 

240.79 

.3039E-4 

5.00 

36.600 

246.56 

.2120E-4 

4.00 

38.240 

250.45 

•1669E-4 

3.00 

40.400 

255.58 

•1227E-4 

2.00 

43.530 

262.96 

-7950E-5 

1.50 

45.790 

268.31 

•5844E-5 

1.00 

49.070 

272.08 

•3842E-5 

.50 

54.650 

265.65 

•1967E-5 

.20 

61.700 

248.11 

•8426E-6 

.10 

66.670 

230.02 

•4544E-6 

I 


I 
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the  significance  of  cloud  position  was  investigated  by  varying  the 
base  of  a  2  Km  thick  cloud  employing  the  drop  size  distribution  of 
the  Deirmendjian  L-Model  cloud. 

For  spherical  particles  the  droplet  volume  per  unit  volumn  of 
atmosphere  is  given  by: 

V  =  ir  /  r^  N(r)dr  .  (5-1 ) 

15  0 

Therefore  the  total  liquid  water  content  per  unit  volumn  of  atmosphere 
is : 

TT  =  npf  r3  N(r)dr  gm  cm"3  (5-2) 

J  0 

_4 

where  r  is  in  units  of  cm,  N(r)  is  in  terms  of  cm  and  the  density  of 
liquid  water  is  1  gm  cm  .  For  the  Deirmendjian  L-Model 

W  =  1 .16678(1 0-2)  gm  cm"2  Km"1  (5-3) 

while  for  the  Marshall -Palmer  size  distribution 

W  =  8 . 8941 5 ( 1 0‘3 )  R'84  gm  cm"2  Km"1  .  (5-4) 

Table  19  contains  the  total  mass  of  liquid  water  within  a  column  of 
2 

1  cm  cross  section  for  the  cloud  models  used. 

5.2  Dependence  on  Layer  Thickness  and  Rainfall  Rate 

To  determine  the  effect  of  varying  cloud  thicknesses  on  the  bright¬ 
ness  temperatures  reaching  the  top  of  the  atmosphere,  cloud  models  were 
inserted  into  the  atmosphere  with  a  constant  cloud  base  of  1  Km.  Thick¬ 
ness  cases  of  1,  2,  3,  4  and  5  Km  were  examined.  Variation  of  bright¬ 
ness  temperatures  over  land  with  respect  to  cloud  thickness  and  rainfall 


-2  -2 

Table  19.  Total  Liquid  Water  Content  (10  gm  cm  ) 


Thickness  (Km)  1  2  3  45 


Dei rmendjian: 
Marshall -Palmer 
R(mm/hr) 

1.16678 

2.33356 

3.50034 

4.66712 

5.83390 

1 

0.88941 

1.77882 

2.66823 

3.55764 

4.44705 

2 

1.59209 

3.18418 

4.77627 

6.36835 

7.96044 

3 

2.23812 

4.47625 

6.71437 

8.95249 

11.1906 

4 

2.84992 

5.69984 

8.54975 

11.3997 

14.2496 

5 

3.43745 

6.87491 

10.3124 

13.7498 

17.1873 

10 

6.15321 

12.3064 

18.4596 

24.6129 

30.7661 

15 

8.65005 

17.3001 

25.9502 

34.0020 

43.2503 

20 

11.0146 

22.0291 

33.0437 

44.0583 

55.0728 

25 

13.2853 

26.5706 

39.8559 

53.1412 

66.4266 

30 

15.4840 

30.9681 

46.5421 

61.9361 

77.4201 

rate  for  channels  1  through  4  are  displayed  in  Figures  3  through  7. 
Figures  8  through  11  depict  the  variation  of  brightness  temperatures 
over  ocean  with  respect  to  cloud  thickness  and  liquid  water  content 
for  the  same  channels.  The  energy  sensed  by  channels  5,  6  and  7  orig¬ 
inates  sufficiently  high  in  the  atmosphere  that  the  brightness  tempera 
tures  are  insignificantly  effected  by  the  cloud  models. 

The  differences  due  to  thickness  for-  channels  3  and  4  may  be 
misleading.  It  is  more  likely  that  a  positional  dependence  for  the 
cloud  is  being  exhibited.  Since  the  cloud  base  was  fixed  at  1  Km, 
only  the  thicker  clouds  reached  into  the  energy  source  regions  for 
channels  3  and  4  (weighting  functions  peak  at  about  7  and  10.5  Km, 
respectively).  For  example,  a  5  Km  thick  cloud  with  a  1  Km  base  at  a 

2  mm/hr  rainfall  rate  results  in  brightness  temperatures  for  channels 

3  and  4  of  237.3  °K  and  228.9  °K.  The  Deirmendjian  L-  Model  cloud 
used  in  the  position  study  has  liquid  water  content  approximately 
equal  to  that  of  the  2  mm/hr  rain  model  (see  Table  19)  and  yet  a  2  Km 


thick  L-Model  cloud  with  base  at  4  Km  results  in  nearly  the  same 
brightness  temperatures  for  channels  3  and  4  as  the  5  Km  thick  cloud 
discussed  above  (237. 2°K  and  228. 7°K,  respectively). 

5.2.1  Over  Land 

Channels  1  and  2  (peaking  near  the  surface  and  2  Km)  show  decreas¬ 
ing  brightness  temperatures  for  increasing  cloud  thicknesses.  The 
decrease  is  greater  for  clouds  modeling  higher  rainfall  rates.  For 
channel  1  at  a  rainfall  rate  of  1  mm/hr  the  decrease  ranges  from  2°K 
for  a  1  Km  thick  cloud  to  16°K  for  a  5  Km  thick  cloud.  For  a  rainfall 
rate  of  30  mm/hr  the  channel  1  brightness  temperature  decreases  14°K 
for  a  1  Km  thick  cloud  and  43°K  for  a  5  Km  thick  cloud.  Channel  2 
displays  the  same  trends  as  channel  1  but  to  a  lesser  degree.  For  a 
rainfall  rate  of  1  mm/hr  the  decrease  in  channel  2  brightness  as  a 
function  of  cloud  thickness  ranges  from  less  than  .5°K  to  5°K.  For  a 
rainfall  rate  of  30  mm/hr  the  range  is  4°K  to  22°K. 

Channels  3  and  4  result  in  relatively  constant  (uneffected)  bright¬ 
ness  temperatures  for  light  rainfall  rates  (less  than  3mm/hr).  For 
rainfall  rates  greater  than  3  mm/hr  the  brightness  temperature  decreases 
with  increasing  cloud  thickness.  The  maximum  decreases  are  7°K  for 
channel  3  and  2°K  for  channel  4;  both  at  rainfall  rates  of  30  m/hr 
and  cloud  thicknesses  of  5  Km. 

5.2.2  Over  Ocean 

The  results  over  the  ocean  are  quite  different  from  the  results 
over  land.  The  only  physical  difference  between  these  two  cases  is 
that  the  surface  emissivity  for  all  channels  is  taken  to  be  .97  for  a 
land  surface  and  .51  for  an  ocean  surface.  This  effects  two  of  the 
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1)  the  surface  contribution  over  ocean  is  approximately  one  half 
of  the  surface  contribution  over  land  (e=.97,  e=.51),  and 

2)  the  atmospheric  contribution  reflected  from  the  earth's  surface 
is  approximately  16  times  greater  over  the  ocean  (r=.03,  r=.49). 

Figures  8  and  9  show  that  channels  1  and  2  follow  similar  trends 

over  the  ocean  but  to  different  degrees.  For  all  cloud  thicknesses  an 

initial  reduction  in  brightness  temperature  is  observed  for  clouds  with 

low  liquid  water  content.  For  channel  1  the  maximum  reduction  is  2°K 

while  for  channel  2  the  maximum  reduction  is  nearly  13°K.  At  some  li- 

*2  -2 

quid  water  content  less  than  5(10  )gm  cm  the  brightness  begins  to 

increase.  The  turning  point  occurs  at  lower  liquid  water  content  for 

thinner  clouds.  The  brightness  then  increases  until  the  liquid  water 
-2  -2 

content  nears  15(10  )  gm  cm  .  For  channel  1  this  increase  results 
in  a  maximum  brightness  over  35°K  higher  than  the  clear  column  value. 
For  channel  2  the  maximum  brightness  is  about  3°K  greater  than  the 
clear  column  value.  For  liquid  water  contents  greater  than  15(10  ) 

O 

gm  cm  brightness  decreases  for  both  channels  and  all  cloud  thick¬ 
nesses.  Note,  the  maximum  brightness  observed  over  the  ocean  is  very 
nearly  the  same  as  the  minimum  brightness  observed  over  the  land. 

Figures  10  and  11  indicate  that  channels  3  and  4  are  relatively 
uneffected  by  a  1  Km  thick  cloud  with  a  cloud  base  at  1  Km.  Increasing 
the  cloud  thickness  above  1  Km  causes  steadily  decreasing  brightness 
temperatures .  For  a  constant  cloud  thickness,  increasing  the  liquid 
water  content  causes  a  decrease  in  brightness.  Note  that  channel  4  is 
effected  by  less  than  1°K  for  all  of  the  cases  studied. 
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5.3  Dependence  on  Layer  Location 

Figures  12  through  23  display  the  importance  of  cloud  position 
relative  to  the  peak  of  the  weighting  function  for  each  channel. 

Channels  1  and  2  show  the  effects  on  a  channel  when  the  cloud  is 
at  or  above  the  peak  of  the  weighting  function,  i.e.,  the  prime  energy 
source  region  for  the  channel.  The  brightness  temperature  decreases 
rapidly  as  the  cloud  is  moved  higher  putting  more  and  more  of  the 
energy  source  for  the  channel  beneath  it.  Eventually,  the  brightness 
becomes  near  constant  as  all  of  the  significant  source  region  for  chan¬ 
nel  energy  is  below  the  cloud  and  therefore  raising  the  cloud  higher 
has  only  slight  effects. 

Channels  3  and  4  show  the  results  of  moving  a  cloud  from  below  the 
energy  source  region  up  through  the  source  region.  Far  below  the 
energy  source  the  brightness  temperature  remains  nearly  constant.  Then, 
as  the  cloud  moves  into  the  source  region,  the  brightness  temperature 
decreases  significantly,  becoming  near  constant  above  the  source  region 
as  for  channels  1  and  2. 

Channels  6  and  7  indicate  that  for  channels  peaking  high  enough  in 
the  atmosphere  to  be  free  of  surface  effects,  brightness  temperatures 
increase  slightly  as  the  cloud  approaches  the  energy  source  region  from 
below. 

Figures  18  through  23  indicate  that  the  same  positional  dependence 
exists  over  the  ocean  as  over  the  land  with  no  change  in  trends  for  any 
of  the  channels. 

5.4  Dependence  on  Profile 

Figures  24  through  31  display  the  results  when  the  sensitivity 
analysis  is  duplicated  for  the  30°N.  Latitude,  July,  profile.  Although 
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the  different  profile  naturally  results  in  different  brightness  tem¬ 
peratures,  the  trends  noted  for  the  N.  Hemispheric  Mid-Latitude  Spring/ 
Fall  profile  are  observed  to  persist. 

5.5  Interpretation  of  Results 

In  order  to  aid  in  the  interpretation  of  results  Figure  32  is  in¬ 
cluded  which  graphically  depicts  the  terms  contributing  to  the  cloud 
boundary  conditions  and  the  mechanisms  for  energy  loss  and  gain  within 
the  cloud.  There  are  three  mechanisms  for  energy  gain  within  the  cloud 
(gaseous  emission,  droplet  emission  and  multiple  scattering)  and  there 
are  three  mechanisms  for  energy  loss  (gaseous  absorption,  droplet  ab¬ 
sorption  and  single  scattering).  The  balance  of  these  mechanisms 
determine  whether  the  emergent  energy  at  the  cloud  top  and  bottom  is 
greater  or  less  than  the  respective  boundary  condition  at  the  cloud 
bottom  and  top.  Of  course,  the  boundary  conditions  themselves  are  an 
important  contribution  to  the  emergent  energy  at  the  opposite  side  of 
the  cloud.  The  upper  boundary  condition  is  totally  uneffected  by  the 
earth's  surface.  For  the  lower  boundary  condition,  quite  the.  opposite 
is  true.  When  the  emissivity  (and  hence  the  reflectivity)  of  the 
earth's  surface  varies,  three  of  the  four  contributing  terms  for  the 
lower  boundary  condition  are  effected: 

1)  the  surface  emission  term, 

2)  the  atmospheric  contribution  reflected  from  the  surface  of  the 
earth,  and 

3)  the  emergent  energy  from  the  cloud  bottom  reflected  from  the 
earth's  surface.  Finally,  it  is  important  to  note  that  the  emission 
terms  are  generally  stronger  near  the  cloud  bottom  due  to  the  tempera¬ 
ture  gradient  within  the  cloud;  and  throughout  the  cloud,  emission  by 
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water  droplets  is  greater  than  emission  by  atmospheric  molecules. 

Introducing  a  cloud  into  the  atmosphere  will  then  result  in  in¬ 
creased  emission  within  the  cloud  layer  and  decreased  transmissivity 
through  the  cloud  layer.  The  three  contributions  to  cloud  top  bright¬ 
ness  are  emission  by  the  earth's  surface,  atmospheric  emission  and 
emission  from  within  the  cloud  layer. 

To  turn  off  the  surface  term,  atmospheric  contribution  or  cloud 
emission  simply  set  the  surface  temperature,  boundary  conditions  or 
cloud  temperature  equal  to  zero,  respectively.  Each  contribution  can 
be  analyzed  separately  by  turning  off  the  other  two. 

Figures  33  through  36  display  the  effects  over  the  ocean  of  cloud 
thickness  and  liquid  water  content  on  the  three  contributions  to  cloud 
top  brightness.  The  solid  line  is  the  graphical  sum  of  these  compon¬ 
ents  and  the  X's  are  computed  cloud  top  temperatures  for  randomly  se¬ 
lected  cases.  Clearly,  the  addition  to  channel  1  energy  by  droplet 

emission  is  greater  than  the  energy  loss  due  to  droplet  absorption  and 

_2 

single  scattering  for  liquid  water  contents  greater  than  .01  gm  cm  . 

This  explains  the  region  of  increasing  brightness  temperature  for 

channel  1  observed  in  Figure  8.  For  liquid  water  content  greater  than 
_2 

.20  gm  cm  the  droplet  emission  reaching  the  cloud  top  is  nearly  con¬ 
stant  and  represents  almost  all  of  the  energy  reaching  the  cloud  top. 
Although  the  surface  emission  and  atmospheric  contributions  to  cloud 
top  brightness  are  small  in  this  region,  they  are  decreasing  with  in¬ 
creasing  liquid  water  content  so  the  cloud  top  brightness  decreases 
slightly  also.  This  explains  the  slight  decrease  in  channel  1  bright¬ 
ness  temperature  in  this  region.  The  fact  that  cloud  emission  accounts 
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for  almost  all  of  the  cloud  top  brightness  for  water  content  greater 
-2 

than  .20  gm  cm  explains  why  the  brightness  temperatures  over  the 
ocean  and  over  the  land  are  nearly  the  same  in  this  region.  The  cloud 
emission  must  be  less  than  the  reduction  in  the  surface  and  atmospheric 
terms  for  liquid  water  content  less  than  .01  gm  cm  in  order  to  explain 
the  slight  decrease  in  brightness  temperatures  in  this  region. 

Over  the  land  the  cloud  emission  curves  remain  the  same.  As  pre¬ 
viously  noted,  the  surface  term  is  nearly  twice  as  large  over  land 
while  the  reflected  terms  are  only  one  sixteenth  as  large.  The  direct 
atmospheric  contribution  remains  the  same.  The  net  result  is  that  over 
land  the  upwelling  brightness  at  the  cloud  bottom  is  significantly 
larger  than  over  the  ocean.  For  example,  the  lower  boundary  condition 
over  land  for  a  1  Km  thick  cloud  ranges  from  280°K  to  285°K  for  rainfall 
rates  of  1  to  30  mm/hr.  Over  the  ocean  the  range  is  200°K  to  240°K  for 
the  same  cloud  thickness  and  rainfall  rates.  Since  the  loss  of  energy 
within  the  cloud  layer  is  proportional  to  the  energy  incident  at  the 
cloud  base,  a  greater  reduction  takes  place  over  the  land.  This  in¬ 
creased  reduction  is  larger  than  the  increase  by  cloud  emission  which 
explains  why  brightness  temperatures  for  channel  1  continuously  de¬ 
crease  with  increasing  liquid  water  content  over  land  (see  Figures  3 
through  7). 

For  channel  2  (see  figures  35  and  36)  the  surface  contribution  to 
cloud  top  brightness  over  the  ocean  is  less  than  that  for  channel  1  due 
to  the  differences  in  transmissivity.  The  atmospheric  terms  are  in¬ 
creased,  however,  due  to  the  shape  of  the  weighting  function  for  this 
channel  and  the  high  reflectivity  of  the  ocean  surface.  The  net  result 
is  a  higher  cloud  top  brightness  for  channel  2.  The  lower  boundary 


condition  is  also  higher  for  the  same  reasons,  resulting  in  a  greater 

reduction  in  the  surface  and  atmospheric  terms  as  they  traverse  the 

cloud.  This  explains  the  deepening  of  the  reduction  in  brightness  for 

n 

liquid  water  contents  less  than  .05  gm  cm  (see  Figure  9).  It  also 

explains  why  the  increase  in  brightness  for  liquid  water  contents  of 

_2 

.05  through  .15  gm  cm  is  not  as  large  as  the  corresponding  increases 
for  channel  1.  Note  that  energy  emitted  by  the  cloud  is  approximately 
the  same  for  both  channels. 

Nearly  the  same  discussion  applies  for  channel  2  over  land  as  for 
channel  1.  The  surface  term  increases  due  to  the  change  in  emissivity 
of  the  earth's  surface.  The  increase  is  smaller  than  for  channel  1, 
however,  due  to  transmissivity  differences.  The  atmospheric  terms  de¬ 
crease  due  to  the  change  in  reflectivity  of  the  earth's  surface.  The 
decrease  is  larger  than  for  channel  1  due  to  the  shape  of  the  weighting 
functions.  The  net  result  is  that  the  lower  boundary  condition  in¬ 
creases  for  channel  2  over  land,  but  the  increase  is  much  less  than  for 
channel  1.  For  example,  for  a  1  Km  thick  cloud  and  rainfall  rates  of 
1  and  30  mm/hr  the  lower  boundary  condition  for  channel  2  is  nearly 
constant  at  285°K  over  land  but  ranges  from  250°K  to  265°K  over  the 
ocean.  Once  again,  this  increase  in  the  lower  boundary  condition  re¬ 
sults  in  increased  reduction  through  the  cloud  resulting  in  continuous¬ 
ly  decreasing  brightness  temperatures  with  increasing  rainfall  rates. 

Channels  3  and  4  are  relatively  uneffected  by  the  emissivity  of 
the  earth's  surface.  This  is  due  to  the  fact  that  the  transmissivity 
from  the  earth's  surface  to  the  top  of  the  atmosphere  is  only  .023  for 
channel  3  and  .003  for  channel  4.  This  explains  why  the  brightness 
temperatures  for  channels  3  and  4  do  not  vary  much  between  land  and 
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ocean  surfaces.  The  highest  cloud  top  modeled  was  at  6  Km.  This  is 
just  below  the  peak  of  the  weighting  function  for  channel  3  (7  Km)  but 
still  well  below  the  peak  for  channel  4  (10.5Km).  As  the  cloud  top 
moves  into  the  energy  source  regions  for  these  channels  emission  by 
the  cloud  is  lower  since  the  temperatures  are  lower  at  higher  altitudes. 
Furthermore,  the  transmissivity  from  6  Km  to  the  top  of  the  atmosphere 
is  less  than  .3  for  both  channels.  Therefore,  within  the  energy  source 
region  occupied  by  the  upper  portion  of  the  cloud,  the  extinction  mec¬ 
hanisms  dominate  over  the  mechanisms  for  adding  energy,  resulting  in 
slight  decreases  in  brightness  temperatures. 

The  interpretation  of  cloud  position  dependence  illustrated  by 
Figures  12  through  23  is  quite  straight  forward.  Note  that  the  lowest 
cloud  base  examined  is  at  1  Km  where  the  brightness  temperatures  over 
ocean  for  channels  1  and  2  are  265. 5°K  and  255. 2°K,  respectively. 

These  values  are  higher  than  clear  column  values  due  to  the  fact  that 
cloud  emission  exceeds  extinction  for  a  2  Km  thick  cloud  with  a  base 
at  1  Km  for  liquid  water  contents  greater  than  .01  gm  cm  (see  Figures 
33  through  36).  The  liquid  water  content  for  the  2  Km  thick  L-Model 
cloud  is  .23  gm  cm  (see  Table  19).  As  the  cloud  is  moved  higher  in 
the  atmosphere  two  changes  take  place.  The  emission  by  the  cloud  re¬ 
duces  due  to  the  reduction  in  cloud  temperature  and  the  lower  boundary 
condition  increases  because  a  greater  portion  of  the  energy  source  re¬ 
gions  for  channels  1  and  2  is  below  the  cloud  base.  This  results  in 
increased  energy  reduction  within  the  cloud  layer.  The  net  result  is 
that  brightness  temperatures  for  channels  1  and  2  decrease  as  the  cloud 
Is  moved  higher  up  in  the  atmosphere.  Eventually,  the  cloud  is  raised 
high  enough  that  nearly  all  of  the  energy  source  regions  for  channels 
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1  and  2  are  already  below  the  cloud  and  cloud  temperatures  are  suffi¬ 
ciently  low  that  emission  by  the  cloud  is  very  small.  Then  moving  the 
cloud  still  higher  has  very  little  effect  and  brightness  temperatures 
are  nearly  constant.  The  same  discussion  applies  over  land  except  that 
over  land  even  the  1  Km  cloud  base  results  in  a  reduction  in  brightness 
temperatures  from  the  clear  column  values. 

Channels  3,  4,  6  and  7  have  weighting  functions  which  peak  higher 
in  the  atmosphere,  reducing  the  effects  of  surface  emissivity.  There¬ 
fore,  the  effects  of  a  2  Km  thick  cloud  are  basically  the  same  over 
ocean  and  land  surfaces.  For  clouds  well  below  the  peak  of  the  weight¬ 
ing  functions  the  cloud  has  no  effect  on  channel  brightness.  This  is 
because  there  is  no  significant  energy  source  below  the  cloud  to  be 
reduced  by  the  cloud  extinction  and  the  transmissivity  from  the  cloud 
top  to  the  top  of  the  atmosphere  is  sufficiently  small  that  emission  by 
the  cloud  is  negligible.  As  the  cloud  moves  into  the  energy  source  re¬ 
gion  a  delicate  trade  off  of  cloud  emission  and  extinction  takes  place. 
For  channels  3  and  4  this  trade  off  holds  the  brightness  temperatures 
very  nearly  constant  until  the  cloud  is  well  into  the  energy  source  re¬ 
gion.  For  channels  6  and  7  a  very  slight  increase  in  brightness  temper¬ 
atures  is  noted.  This  increase  is  less  than  .2°K,  however,  which  is 
below  the  noise  level  of  the  SSM/T.  Finally,  as  the  cloud  is  moved 
above  the  peaks  of  the  weighting  functions  brightness  temperatures  de¬ 
crease,  eventually  becoming  near  constant  as  for  channels  1  and  2. 
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Figure  5.  Brightness  temperature  vs  rainfall  rate  Figure  6.  Brightness  temperature  vs  rainfall  rate 
(or  liquid  water  content)  for  3  Km  thick  (or  liquid  water  content)  for  4  Km  thick 

cloud  over  land  (Mid-Latitude  Spring/  cloud  over  land  (Mid-Latitude  Spring/ 

Fall  profile).  Fall  profile). 


9  Km 


Brightness  temperature  vs  rainfall  rate  (or  liquid  water  content)  for  5  Km  tnic'  cloud  over 
land  (Mid-Latitude  Spring/Fall  profile). 
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Figure  8.  50.50  GHz  brightness  temperature  as  a  function  of  cloud 

thickness  and  liquid  water  content  (Mid-Latitude  Spring/ 
Fall  profile  over  ocean). 


Figure  9.  53.20  GHz  brightness  temperature  as  a  function  of  cloud 

thickness  and  liquid  water  content  (Mid-Latitude  Spring/ 
Fall  profile  over  ocean). 
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Figure  10.  54.35  GHz  brightness  temperature  as  a  function  of  cloud 

thickness  and  liquid  water  content  (Mid-Latitude  Spring/ 
Fall  profile  over  ocean). 


Figure  11.  54.90  GHz  brightness  temperature  as  a  function  of  cloud 

thickness  and  liquid  water  content  (Mid-Latitude  Spring/ 
Fall  profile  over  ocean). 
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Figure  16.  58.82  GHz  brightness  temperature  Figure  17.  59.40  GHz  brightness  temperature 

dependence  on  cloud  base  height  for  dependence  on  cloud  base  height  for 

a  2  Km  thick  Deirmendjian  L-Model  a  2  Km  thick  Deirmendjian  L-Model 

cloud  (Mid-Latitude  Spring/Fall  cloud  (Mid-Latitude  Spring/Fall 

profile  over  land).  profile  over  land). 


50.50  GHz  brightness  temperature  Figure  19.  53.20  GHz  brightness  temperature 

dependence  on  cloud  base  height  for  dependence  on  cloud  base  height  for 

a  2  Km  thick  Deirmendjian  L-Model  a  2  Km  thick  Deirmendjian  L-Model 

cloud  (Mid-Latitude  Spring/Fall  cloud  (Mid-Latitude  Spring/Fall 

profile  over  ocean).  profile  over  ocean). 


Figure  20.  54.35  GHz  brightness  temperature  Figure  21.  54.90  Ghz  brightness  temperature 

dependence  on  cloud  base  height  for  dependence  on  cloud  base  height  for 

a  2  Km  thick  Deirmendjian  L-Model  a  2  Km  thick  Deirmendjian  L-Model 

cloud  (Mid-Latitude  Spring/Fall  cloud  (Mid-Latitude  Spring/Fall 

profile  over  ocean).  profile  over  ocean). 
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50  GHz  channel  profile  dependence  over  ocean. 


channel  profile  dependence  over  ocean 
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Figure  32.  Radiative  transfer  through  a  cloud  layer 


Figure  35.  Component  contributions  to  Figure  36.  Component  contributions  to 

53.20  GHz  cloud  top  brightness  53.20  GHz  cloud  top  brightness 

for  a  1  Km  thick  cloud  over  for  a  3  Km  thick  cloud  over 

ocean  (Mid-Latitude  Spring/  ocean  (Mid-Latitude  Spring/ 

Fall  profile).  Fall  profile). 


CHAPTER  6 


TEMPERATURE  RETRIEVAL 

The  difficulty  in  reconstructing  a  temperature  profile  from  ra¬ 
diance  or  brightness  temperature  values  is  due  to  the  fact  that  the 
Fredholm  equation  with  fixed  limits  may  not  always  have  a  solution  for 
an  arbitrary  function.  The  extraction  of  T(z)  from  under  the  integra¬ 
tion  sign  of  Equation  (4-19)  is  possible  because  for  certain  frequency 
ranges  the  integrand  reaches  a  strong  maximum  at  different  heights. 

6.1  AFGWC  Statistical  Method 

The  SSM/T  retrieval  method  is  based  on  a  linear  multiple  regres¬ 
sion  technique.  It  is  assumed  that  the  deviation  of  the  parameter 
vector  (temperature  profile)  from  the  climatological  mean  may  be  ex¬ 
pressed  as  a  linear  combination  of  the  deviations  of  measured  data 
(brightness  temperatures)  from  its  mean.  Thus: 

p  -  <p>  =  D(d  -  <d>)  (6-1) 

A 

where  p  is  the  atmospheric  temperature  vector,  p  is  the  predicted 
value  of  p,  <p>  is  the  climatological  mean  of  p,  d  is  the  data  vector 
of  brightness  temperatures,  <d>  is  the  vector  of  brightness  temperature 
mean  values,  and  D  is  a  matrix.  The  0-matrix  is  determined  by  the 
condition  that  the  square  of  the  difference  between  the  true  value  of 
the  vector  p  and  the  predicted  value,  p,  be  a  minimum.  This  leads  to 
the  equation: 
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D  =  C(p-<p>,  d-<d>)  C'1(d-<d>,  d-<d>)  (6-2) 

where  C(.,.)  denotes  the  covariance  matrix  of  two  arguments.  Eq.  (6-1) 
can  be  written  as: 

p  =  Dd  +  A  (6-3a) 

where 

A  =  <p>  -  D<d>  (6-3b) 

so  that  the  predicted  value  of  the  temperature  vector  is  expressed  as 
a  constant.  A,  and  a  linear  combination  of  the  data. 

The  covariance  matrix,  C,  required  in  Eq.  (6-2)  has  been  deter¬ 
mined  from  apriori  data  for  impl imentation  purposes.  That  is,  the 
microwave  response  for  a  large  number  of  archived  atmospheres  repre¬ 
senting  a  wide  range  of  meteorological  conditions  was  calculated,  and 
from  this  information  the  D-matrices  were  determined.  After  impl iment¬ 
ation,  real  time  data  will  be  used  to  update  these  D-matrices  on  a  re¬ 
gular  bases. 

Since  microwave  radiometers  posses  inherent  instrumental  noise, 
this  noise  must  be  accounted  for  in  the  computations.  The  data  vector, 
d,  may  be  written  as: 

d  =  dQ  +  s  (6-4) 

where  dQ  is  the  exact  brightness  value  and  6  is  a  random  error  repre¬ 
senting  system  noise.  Then, 

C(p-<p>,  d-<d>)  *  C(p-<p>,  dQ-<d0>)  ( 6-5 ) 

and 


4 
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C(d-<d>,  d-<d>)  *  C(d0-<d0>,  dQ-<d0>)  +  C (6,6)  (6-6) 

Eq.  (6-6)  incorporates  the  instrumental  noise  characteristics  in  the 
computations  through  the  noise  covariance  matrix,  C(<5,5). 

Eq.  (4-29)  may  be  rewritten  to  represent  clear  column  brightness 
at  the  top  of  the  atmosphere  as: 

oo 

Tgfv,*^)  *  £v^sTv(0»“'*vji)  +  /  T(z)dTv(z,»;u.j) 


+  (1-e  )t  (0,»;y.)[/  T(z)dx  (0,z;-u.)]. 

v  V  *0  ^ 


(6-7) 


Rewriting  Eq.  (6-7)  to  obtain  an  expression  for  e^Ts  yields: 


evTs  =  CTb(v ,-iWi >-xv(0.-iw1 ) 0 -cv)  /  T(z)dTv(0,z;-Ui) 


00 


-  /  T(2)dT  (z,®iy.)]/T  (0,®;y.) 

0  V  1  V  i 


(6-8) 


Now  defining: 

Tsky  =  l  T(z)dTv<0.z;-Ul)  , 

00 

Tatn  =  l  T(i)dTv(z.-i,,)  . 

and 

fg(^)  -  Tg('J»"*v^)  » 

Eq.  (6-8)  may  be  expressed  more  simply  as: 

eJs  =  CTB(',)‘Tv(0’a,;wi)(1'ev)T  L  “T  +  3/t  (O.cojp.)  . 

sky  atm 


* 


Then 
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e  T  t  (0,‘ 
v  S  v 


’Jy^  )(1 


~f^  =  tTB(v^‘Tv(0, 


>,Mi ^Tsky 


Tatm^ 


and 


e  T 

v  S 


*  [T8tv)-Tv(0.-iu1)T$ky-T1|ttl|]/CTv(0..;.1)O-  -p)]  (6-9) 


In  particular. 


(6-10) 


Eq.  (6-10)  is  used  to  obtain  an  expression  for  TyTs  in  terms  of  the 
channel  1  brightness. 

For  DMSP  application  a  data  vector  was  chosen  to  be  dependent  only 
on  atmospheric  conditions.  That  is: 


d.  =  t  (0,°°)  T  i  +T  . 
l  vv  '  sky  atm 


From  Eq.  (6-9)  we  may  write: 


(6-11) 


(0’"i'‘1>WTatn,  -  -f?>  =  d1  (6‘12) 


Since  the  terms  on  the  left  are  clearly  atmospheric  dependent  only, 
the  correction  to  brightness  temperature  represented  by  the  negative 
term  on  the  right  will  clearly  remove  the  effects  of  surface  emission 
and  reflection  from  the  brightness  value.  Substituting  Eq.  (6-10) 
into  Eq.  (6-12)  yields: 


di  -  Tb(v)  -CTb(50.5)-S50  5]aVi 


(1*1,  2, • • • ,  6) 


(6-13) 
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where 


^50.5  "  ^Tv^’“,lJi ^sky  +  ^atm-^g  ^ 


and 


[Tv(o..i,1)o-  -^nVj 

lv.  =  Y~ 

[Tv(0,«;Wi)(1- -^)]50> 


The  values  for  s^q  g  and  the  ay  have  been  computed  for  implimentation 
from  the  same  apriori  data  used  to  compute  the  D-matrices.  Then  the 
multiple  regression  Equations  (6-3a)  and  (6-3b)  may  be  expressed  as: 


P  ’  STb  +  A’ 
and 


(6-1 4a ) 


A'  =  <p>  -  D<d>  +  s^q  5  D  a  (6-1 4b) 

where 


and  S=(-Da:D)  is  a  matrix  whose  first  column  is  -Da  and  whose  remain¬ 
ing  columns  are  the  columns  of  D. 

6.2  Chahine's  Relaxation  Method 

As  a  consequence  of  the  mean  value  theorem,  the  contribution  to 
brightness  temperature  due  to  atmospheric  emission  can  be  expressed  as 
the  change  in  upwelling  brightness  at  one  level,  z.  (the  level  where  the 

J 

channel  weighting  function  peaks),  times  the  effective  width  of  the 


<* 


1 
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weighting  function,  a-z.  Then 

J 


tb(v°  ‘  VvTvt0,*>  *  T(zj)[-~!rL\  y 


(6-15) 


where  Ts  and  T(z j )  are  the  true  temperature  values  and  Tg(vj,®)  is  the 
observed  brightness  temperature.  If  we  arbitrarily  guess  a  temperature 
profile,  then  the  calculated  brightness  temperature  can  be  represented 
as: 


3t  1  (z ,®) 


-  T'(zj)t-  3Z 

Dividing  Eq.  (6-15)  by  Eq.  (6-16)  yields: 


-] 


4iz 


TB(vj,-)-TscvTv(°.») 


3t  (z,®) 

_  J 

3t '  (z ,®) 

]  ,  y 

j 


(6-16) 


(6-17) 


Since  the  transmissivities  and  weighting  function  values  do  not  change 
very  much  as  a  function  of  temperature. 


(z,®) 

v 

I- ajz 

J 


=  1 


Therefore, 


T‘(z)[TB(v1,®)-TsVv(0,®)]  a 

CTB  (vj»A)‘TssvTv(°»“)] 


(6-18) 


When  the  contribution  to  brightness  temperature  at  the  top  of  the 
atmosphere  due  to  the  surface  emission  term  is  dominate  or  negligible. 


we  obtain  the  relaxation  equation  of  Chahine: 


T‘  ( 2  )Tg  ( v  j  >°°) 

TB(vj*”) 


T(z) 


(6-19) 


The  itterative  method  of  solving  for  the  correct  temperature 
values  corresponding  to  the  peak  levels  for  the  weighting  functions  of 
j  channels  is  as  follows: 

1)  Make  and  initial  guess  (n=0)  for  the  Tn(z.),  where  the  super- 

J 

script  denotes  itteration  number. 

2)  Interpolate  temperature  values,  Tn(z^),  to  N  standard  levels 
from  the  Tn(z.). 

J 

3)  Using  an  accurate  quadrature  formula,  evaluate  the  correspond¬ 
ing  brightness  temperature  values,  Tg(vj»°°)  using  Eq.  (4-17). 

4)  Compare  the  computed  brightness  temperatures  with  the  observed 
values,  TD(v . ,») .  If  the  residuals, 

D  j 

Rj  =  |TB(vj,»)-TB(vj,») |/IB(vj,»)  , 

are  less  than  a(a-*0)  for  each  channel,  then  Tn(zk)  is  a  solution. 

5)  If  any  of  the  residuals  are  greater  than  a,  then  Eq.  (6-19)  is 
applied  to  each  channel  to  obtain  a  new  estimate  of  the  temperature  for 
the  level  corresponding  to  the  peak  of  the  weighting  function  for  each 
channel . 

6)  Go  back  to  step  2)  and  repeat  until  all  residuals  are  less 
than  a. 


6.3  Presentation  of  Results 

The  two  retrieval  methods  described  in  Sections  6.1  and  6.2  were 
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used  to  retrieve  the  Mid-Latitude  Spring/Fall  profile  used  in  this 
sentivity  study.  The  results  of  these  retrievals  are  depicted  in  Fig¬ 
ures  37  through  40.  In  each  case  the  curve  labeled  (a)  is  the  actual 
profile  used  to  calculate  the  brightness  temperatures.  The  curve 
labeled  (b)  is  the  retrieval  accomplished  using  clear  column  brightness 
temperature  values.  Curves  labeled  (c)  through  (f)  represent  precipi¬ 
tating  cloud  cases  of  increasing  rainfall  rates  and  thicknesses. 

Figure  37  and  38  indicate  that  the  AFGWC  method  has  very  success¬ 
fully  removed  surface  effects  from  the  retrieval.  Note  how  similar  the 
retrieved  profiles  are  over  land  and  over  ocean  for  the  same  thickness 
and  rainfall  rate.  Curve  (c)  on  both  figures  represents  a  light  to 
moderate  precipitation  case  (up  to  5  mm/hr  and  2  Km  thickness).  The 
surface  temperature  suffers  some  degradation  from  the  clear  column  case 
but  still  stays  within  10°K  of  truth.  Above  3  Km  the  retrieval  for  all 
levels  is  within  .5°K  of  the  clear  solumn  retrieval.  This  indicates 
that  for  light  to  moderate  rain  and  certainly  for  nonprecipitating 
cloudy  atmospheres,  the  AFGWC  method  will  provide  a  global  data  base  of 
atmospheric  temperature  profiles  representing  a  substantial  improvement 
over  climatology.  The  method  does  perform  significantly  worse  for 
heavier  rainfall  rates  and  thicker  clouds.  The  worst  errors  occur  at 
the  surface  and  near  the  tropopause.  Rainfall  rates  of  this  magnitude 
(greater  than  5  mm/hr)  do  not  occur  over  very  large  areas  or  for  very 
long  time  periods.  Therefore,  they  should  not  have  a  significant 
effect  on  global  mapping,  although  they  do  point  out  the  need  for 
careful  quality  control  of  individual  profiles. 

Figures  39  and  40  indicate  that  Chahine's  relaxation  method  per¬ 
forms  very  well  below  the  tropopause  for  clear  atmospheric  conditions. 


temperature  retrieval  over  ocean  (Mid-Latitude  Spring/Fall 


temperature  retrieval  over  land  (Mid-Latitude 


Over  land  the  surface  temperature  becomes  poor  very  rapidly  for  even 
moderate  precipitating  conditions.  For  a  2  Km  thick  cloud  with  a  pre¬ 
cipitation  rate  of  5  mm/hr  the  retrieved  surface  temperature  is  19°K 
too  low.  The  quality  of  the  retrieved  temperature  profile  continues 
to  degenerate  rapidly  as  the  degree  of  precipitation  increases. 

Over  the  ocean  Chatline' s  method  continues  to  do  well  for  clear 
atmospheric  conditions.  At  a  rainfall  rate  of  1  mm/hr  with  a  1  Km 
thick  cloud,  the  temperature  retrieval  begins  to  worsen  in  the  lower 
troposphere,  although  the  surface  temperature  retrieved  is  good.  For 
a  2  Km  thick  cloud  and  a  5  mm/hr  rainfall  rate,  this  method  fails 
completely,  with  the  retrieved  profile  oscillating  wildy  about  the 


truth. 


CHAPTER  7 


CONCLUSION 

Channels  1,  2,  3  and  4  of  the  SSM/T  have  been  shown  to  be  effected 
by  precipitation  to  various  degrees. 

For  channels  3  and  4  the  effects  are  quite  small  and  in  fact  are 
of  an  order  close  to  that  of  the  instrument  noise  level.  In  other 
words,  the  variations  noted  for  these  channels  are  basically  insignifi¬ 
cant  to  operational  applications. 

Channels  1  and  2  display  a  much  more  significant  response  to  pre¬ 
cipitation.  General  trends  in  the  brightness  variations  for  these 
channels  are  dependent  on  cloud  liquid  water  content,  thickness  and 
surface  emissivity. 

For  a  theoretical  retrieval  scheme  based  on  the  premise  that  the 
input  brightness  temperature  values  are  clear  column  values,  these 
variations  cause  a  rapid_ degradation  in  the  accuracy  of  the  retrieved 
temperature  profile.  However,  methods  are  available,  such  as  the  AFGWC 
statistical  method,  which  can  somewhat  adjust  to  surface  and  lower 
tropospheric  conditons.  When  clear  column  conditions  exist,  the  re¬ 
trieval  can  be  made  without  degradation  due  to  surface  emissivity.  For 
precipitating  conditions,  the  surface  emissivity  is  underestimated,  but 
this  error  tends  to  compensate  for  extinction  within  the  precipitation 
layer  reducing  the  brightness  error  in  channels  sensing  the  lower 
troposphere.  Although,  very  high  rainfall  rates  can  result  in  a  very 
poor  temperature  retrieval,  the  frequency  and  spatial  coverage  of  such 
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atmospheric  conditions  are  sufficiently  low  that  global  mapping  of 
vertical  temperature  profiles  derived  from  microwave  sensing  from 
satellite  will  provide  a  significant  improvement  over  climatology 
for  forecast  model  initializations. 


i 
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